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Random tensor models are generalizations of matrix models which also support a 1/N expansion. 
The dominant observables are in correspondence with some trees, namely rooted trees with vertices 
of degree at most D and lines colored by a number i from 1 to D such that no two lines connecting 
a vertex to its descendants have the same color. In this Letter we study by independent methods a 
generating function for these observables. We prove that the number of such trees with exactly pi 
lines of color i is ( E *=i Pi+1 ) ■ ■ ■ ( E ?=i p ' +1 ) . 
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I. INTRODUCTION 



Qs, ' The study of large random matrices in the past thirty years has successfully described measures which can be 
written as the exponential of single trace invariants perturbing a Gaussian. In addition to the standard Fcynman 
diagrammatic expansion [l| , powerful methods exist to solve such models 0] , including orthogonal polynomials (which 
rely on eigenvalue decomposition) and more recently the topological expansion developed by Eynard The latter 
starts with the Schwinger-Dyson equations (also known as loop equations) written in terms of the resolvent Tr z _^ IM - t 
• where M is the random N x N matrix, and provides an intrinsic way to solve them at all orders in the 1/N expansion. 

Random tensors are a generalization of random matrices to rank D objects (having D indices of size N each). The 
study of their large size statistics became possible thanks to the 1/N expansion discovered in Q. This expansion 
relies on the construction of multi-unitary invariants, i.e. tensor contractions which are invariant under the external 
'— — t ' tensor product of D copies of U(N), each of them acting independently on each tensor index. In contrast with random 
matrices, there are many invariant monomials at a fixed degree. These monomials are indexed by colored graphs @, H| 
^ ' (these colors correspond to the position of the index in a tensor contraction: in T ai ... aD1 a\ has position, hence color, 
, 1, up to ar> which has color D). The 1/N expansion generalizes to any measure on a single random tensor that can 
• be written as the exponential of such invariants [f|. 
£S| \ In the diagrammatic approach, the graphs contributing at leading order, known as melonic graphs (equivalent to 
■ planar graphs of matrix models), have been identified in [6[ enabling one to solve these models exactly at leading order 
\Q \ Q. At the core of these solutions, a universality theorem, first derived in Q, states that all models are Gaussian 
. at large N (but with a covariancc which crucially depends on the joint distribution). In particular all invariant 
O^l ' monomials corresponding to melonic graphs at fixed degree have the same expectation value. However, this certainly 
does not hold at sub- leading orders in the 1/N expansion, and different melonic graphs have different sub- leading 
corrections in 1/N. 

A method which would be fruitful to adapt to tensor models is the one developed by Eynard Q . It first requires to 
introduce an equivalent of the resolvent. Writing the matrix resolvent like Tr z _^ IMi = J2n>o z ~ n ~ 1 Tr(MMt)", it is 
tempting to generalize it by changing Tr(MM^)™ with the sum over all invariants of degree n. As z counts the degree 
of each invariant, this object does not distinguish different invariants having the same degree. Although sufficient 
for the study of the leading order, this object should be refined in order to explore the finer structure of sub-leading 
orders, for which a generating function which probes the structure of each invariant beyond their degree seems better 
adapted. 

This Letter is a modest contribution in this direction. The dominant melonic invariants are in one-to-one correspon- 
dence with D-aiy trees with lines colored 1 up to D such that no two lines connecting a vertex to its descendants have 
the same color. At leading order only the number of vertices matters, but at subsequent orders one must distinguish 
between different colored trees having the same number of vertices. As a first step, we study the generating function 
of these colored trees and find an explicit formula counting how many such trees with pi lines of color i one can build. 
Countings of colored trees exist in the literature (like Q) but we could not find the counting which is relevant for our 
purpose. 
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In Section [II] we introduce the problem and state our main results. The generating function is presented in Section 
IIIII and the proof of our results, based on a linear recursive sequence related to the generating function, is contained 
in Section llVl 

II. STATING THE PROBLEM 

We consider a family of rooted trees defined by the following properties 

• each vertex has at most D descendants, D > 2, 

• each line receives a color i = 1, ...,D such that no two lines connecting a vertex to its descendants have the 
same color. 

By adding leafs (univalent vertices) appropriately any such tree becomes a rooted D-ary tree with colored lines. Let 
C Pl ... PD be the number of such trees with exactly pi lines of color i for all z = 1, ... ,13. The purpose of this note is 
to derive an explicit formula for C Pl ... PD . 

Our strategy is based on the generating function 

oo D 

F( 91 ,..., gD )= j2 c Pi- PD n»r. (1) 

Pi,...,Pd=0 i=l 

and the sequence {F n (gi)) n ^ defined by 

F = l, F n { gi ,...,g D )= J2 <%.. w f[g? with C ;\... pD = n f[ ( E ^ Pi + ") . (2) 

P!,...,pn i=l X-,i=lPt~T~ n j=l ^ ' 

We will prove that (F n ) is a linear recursive sequence whose characteristic equation is an algebraic equation satisfied 
by F . Examining the roots of this algebraic equation, we will obtain the following proposition. 

Proposition 1. The sequence {F n {gi)) n ^ is a geometric sequence with common ratio F, 

F n ( gi ,...,g D ) = (>( ffl ,..., . (3) 
This implies that F(gi, . . . ,.gr>) = F\{g\, . . . ,go), hence the following corollary. 

Corollary 1. The number C Pl ... PD of rooted line-colored trees with maximal degree D and exactly pi lines of color 
i = 1 , . . . , D is 



III. THE GENERATING FUNCTION 

The generating function ([1} satisfies an algebraic equation which is obtained by simply observing that the root of 
a tree can have k < D descendants connected by lines of colors z 1; . . . , i). all different. Therefore 



F( gi ,..., 9D ) = l + ( 9h) F +( E 9i 2 9i 2 )F 2 + --- + (gi---g D )F D , 

(5) 



l<ii<D l<n<i 2 <D 
D / k 



E E j9i l )[F(gi,...,g D )} h 



k=0 s l<h<---<i k <D 1=1 

and F is the root 2Vo)(<?i) ■ • - .9d) of this polynomial equation such that 

lim x {0) { gi ,...g D ) = 1 . (6) 

The following lemma show how to distinguish the desired root iE(o)(<7i, . . . go) from the other roots of the above 
polynomial. 
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Lemma 1. For all R > 1, there exists e R > such that for all \gi\, ■ ■ ■ , \9d\ < £r, the polynomial 

D k 

q(x) = -x + y / ( £ U^) xk (7) 

fc=0 l<ii<".<t*<£) (=1 

has exactly one root X(o) with |x(o) | < -R, other roots xn\ for i ^ having norms \xu-\ \ > R. In particular X( ) = -F 1 
as it is the only root satisfying lim gi _>o 2J(o) = 1. 

Proof. To establish this lemma we use Rouche's theorem whose statement is now recalled. Let /, g be two 
holomorphic functions and S a closed contour which does not contain zeros of / and g. If for all z € S 

\f(z)-g(z)\<\g(z)\, (8) 

then the number of zeros of / and the number of zeros of g encircled by S (counted with multiplicities) are the same. 
We now exploit this theorem to get bounds on the absolute value of the roots of the polynomial (J7]). Take 

f(z) = Q(z), and g{z) = -z . (9) 

Let R > 1, and S — {z G (D, \z\ = R} be the circle of radius R. Note that on S, \g(z)\ = R. Set e R > with 
eji < R 1 R ~ x so that on S, 

\f(z) - g(z)\ < £( £ ' ' ' toj) l^| fe < £ (fc) ^ C 1 + e «^) D ^ ■ ( 10 ) 

fc=0 l<i a <i 2 ---<i fc <ZJ fc=0 ^ ' 

As g(z) has an unique root z — inside the circle of radius R, we conclude that f(z) also has exactly one root, which 
we denote xrm, with |:£(o)| < R- As this is the only root which can go to one when gi — > 0, it is identified with the 
generating function F. □ 

Remark 1. At D = 2, Q(X) = 1 + (gi + g2 — 1)X + gig^X 2 and Q{x) = is easily solved, 



1 - .91 - g-2 - - gi - .92) 2 - 4.9152 , 1 - .91 - .92 + VC 1 ~ -9i - .92) 2 - 4gig 2 1 > 

x (0 ) = ^— and = — . (11) 



IV. THE LINEAR RECURSIVE SEQUENCE 



First we show that the sums defining each F n (gi) in ([2]) converge absolutely when < - — jjh for all t = 1, . . . , D. 

D-l 



Let e > such that e < ^ D *)> , then for all complex gi,. . - ,gD with norm |<jj| < e 



p=o { Ps } 

J2i Pi=P 



\F n { 9l ,..., 9D )\<^-^ £ ( P + n )-( P + n ) (12) 
^ p + n j—* \ Pi J V PD 1 



The sums over pi are computed by equating the coefficients of x p in (1 + x) Dp+Dn and (l + x) p+n ■ • • (1 +x) p+ ™, hence 

i^i,..„9.)i<E^r p+ ^)^- ( 13 ) 



One finds using the Stirling formula that the radius of convergence of the above series is - — jjb > e, hence 

F n (gi, . . . , gu) converges absolutely. 

Lemma 2. The sequence (F n ) respects the recursion 

D 



Vn > F n+1 = F n + £ £ 1 F n+k ■ (14) 



k=l x l<h< - <i k <D 1=1 
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Proof: The recursion translates into 

D 

VPi > 1 C PU 1 ..PD = C PU-PD + E E C 'pit A '.p. 1 -l,...p, fc -l,...PB ' ( 15 ) 

l<i 1 <i 2 ---<ifc<.D 

The boundary cases, when some pi = 0, just reproduce the recursion at level D — 1. Let us denote P = y~]? 1 Pi- The 
right hand side of (|T5|) is 



(P + n)"^(P+_^] 
' - Pi + n -+ 



J(P- Pi + n+ l) + J2(n + k) Yl IP*. (P-Pi + n+1) . (16) 



rc=ift!0P -ft + n + l)! 



i=l fc=l l<ii <-<i k <D 1=1 jjti 1: ...,i k 



We write n Oi=i (-f — Pi + n + l) = (n + 1) IL^i^ — P» + n+l) — ]^[ i=1 (P" — pi + n+ 1) so as to re-arrange the square 
bracket above as 



z fc— 1 l<«i <-"<ifc <.D / — 1 j^i±,...,i k 

D D D k 

^(n+^HiP + n+^-HiP-p. + n+^ + ^ik-l) E 11^ II (i'-Pi+n+l) 

Z=l 2=1 fc = l l<ii <■ ■ ■ <Zfc <L> /=1 

= (n + 1)(P + n)(P + n+ l) " 1 + (n + 1)(P + n + l) ' 1 - ]J(P - Pl + n+l) 



+ X>-1) E 11^ II (P-Pj+n + 1). (17) 

k=l l<ii< - <i fc <D ;=1 j^ii,...i k 

The first term of the last equality is exactly what is needed to form C'™ 1 + . 1 pD . Therefore we focus now on the sum of 
the three other contributions, 

D D k 



(n+l^P + n+l) - 1 -H(P + n+l- Pl ) + J2(k-l) ]T JJpi, ]J (P + n+l- Pj ) 

i=l k=2 l<i 1 <-<i k <D 1=1 j^i u ...i k 

D k 

= {n + l){P + n + l) D - 1 -{P + n+l) D + {P + n + l) D - 1 P + Y J {-) k+1 E J^]( P + n + 

k=2 l<ii<---<j fc </5 1=1 

D k D-k m 

+E( fc ~ 1 ) E ftp*. E( p +^+ 1 )^ m (-r E Up*- ( 18 ) 

k=2 l<h< - <i k <D 1=1 m=0 l<ji < - <j" m <D t=l 

The first three terms cancel. For the remaining terms 

ER fe+i e [n^](^+-+i) D - fc 

k=2 l<ii<-<i k <D 1=1 

D k D-k m 

+E( fc -!) E E( p +" +1 )^ m (-r E n^- ( 19 ) 

k=2 l<ii<-<i k <D 1=1 m=0 1<j"i < — <j m <D t=l 

we take into account that k + m ordered integers can be partitioned into ( fc ~£ m ) ways into two subsets of k and m 
ordered integers. Thus the second sum rewrites as a sum over q = k + m 

D q , 

<1 



E E Pn---P^H g (P + n + l) D -' 1 J2(-) k ( k - 1 )(l) ■ ( 2 °) 

q=2 l<ii<-<i,<£) k=2 



But 



EH fc (fc - i) 1 1) = i + EH fc (fc - 1) [lj = i - (i - 1) 9 - 9(1 - I)*" 1 = i , (21) 



k=2 x 7 k=0 
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hence the whole quantity displayed in (|T9|) is zero and the lemma [2] follows. □ 

Therefore, F n+ r> is obtained recursively from the set (F p ) p<n+ D- The characteristic polynomial of this recursion 
is exactly Q(X) (Equation ©). It means that the solution xrfy is one of the common ratios of (F n ). We denote the 
others xu^, and assuming they are all different, 

F n = ax? 0) +J2bj%i (22) 

3 

for some functions a(gi, . . . , go), bj(gi, . . . , go) which can in principle be determined by D initial conditions. However, 
we cannot use the initial conditions (remember we want to prove that F = F\) so we have to proceed differently. 
Each common ratio in the sum (|22j) is controlled thanks to the Lemma [IJ as x^ is bounded from above and each 
xu\ is bounded from below. Now we need to control the sequence (F n ) independently of its common ratios. This is 
done through the following lemma. 

Lemma 3. For all K > 1, there exists ck > such that for all g±, . . . , gr> G C with \gi\, ■ ■ ■ , \9d\ < e K, F n {g\, . . .go) 
is polynomially bounded by K , 

Vn>0 \F n ( 9l ,...g D )\ < K n . (23) 
Proof. For n = 0, this is trivial as Fo = 1. Let thus be n > 1 and K > 1. It is enough to chose ck such that 

(D-l) - 1 1 9nc 1 2Dee K , x 

(D — 

With < * — ^jb , one can use Equation (fT5|) which implies 



p=0 p=0 



(25) 



We use the fact that (n + p) p < (2n) p when p < n and (n + p) p < (2p) p when p > n to obtain the bound 

< J2n p[ -^- + £ ^(2De K r < e 2D " e - + £ P -(2De K f . (26) 

p— p=n p— n 

Now we use p! > ^2iip e plnp_p , Vp > 1 and as ex < we get 

1 oo 1 oo 

\F n \ < e 2DneK + —Y— (2Dee K ) p < e 2Dn * K + ^= Y{2Det K ) p 

^UVp ^ (27) 

< e 2Dne K , 1 2Jee A - < f 2DtK 1 2J>ee Jf \ n 

V2^ 1 - 2Z>ee jr ~ V ^\-2Dee K ) 

□ 

We are now in position to prove Proposition [TJ by combining Lemmas [I] and [3J Choose ii > K > 1 and consider 
\gt\ < inf(ejj, ejf) with eR,ex as in the lemmas. First we consider the case where all X(j) have different norms and 
denote xr max ) ^ the one with the largest norm. In particular |x( max )| > R > 1. At large n, the norm of F n is 
dominated by b/ max \Xf max \. Hence there exist a constant A > and an integer N such that for all n > iV, 

l-FnCffi, ■ ■ • ,5d)| > A |6 (max) | |rr max |™ > A \\ max) \R n . (28) 

However |F n | < K n with i? > K. Therefore we conclude br max \ = 0. We can repeat this reasoning with the root xu\ 
that has the second largest norm, and so on until we get F n = ax^ . The initial condition Fq = 1 for all gi finally 
leads to F n = x? y 

The case where some of the roots have the same norm is quite similar. The idea is to extract sub-sequences (F r ^) 
for which F r r n \ behaves at large n like a coefficient times some combination of the roots j ^ 0, where this 
combination is greater than R n when \xu\ \ > R. 

Remark 2. At D = 2, the number of line-colored trees with p\ lines of color 1 and pi lines of color 2 is C PlP2 = 
1 (P1+P2+ ^pi+P2+ ^ These numbers are known as the Narayana numbers N(p\ + p2 + l,Pi + !)■ 
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Remark 3. By summing the numbers C Pl ... PD over all possible numbers of lines of each color at a fixed total number 
of lines P — 1, 

V C 1 (D(P-1) + D\ 1 (DP + 1\ , . 

pl " PB ~P\ P-l J-DP+l{ P )' [J) 

we obtain the total number of D-ary trees on P vertices also known as the D-Catalan numbers (pp. 200 in fldi] . 
proposition 6.2.2 in F7H / and more details in J&J). 

Remark 4. Proposition^ implies that F n F m = F n+m , corresponding to interesting combinatorial identities, 

{fei=0,..., Pi } 

For example, when D = 2 one gets 

P ^-? n m fk\ + k 2 + n\ fk\ + fc 2 + n 

u T^n + ^2 + n) (pi - h+p 2 - k 2 + m)\ k\ )\ k 2 

Pi — k\ + p 2 — k 2 + m\ [pi — k\ + p 2 — k 2 + m\ n + m fpi+P2 + m + n\ (pi+p 2 +m + n 



p\ — k\ I \ p 2 - k 2 I pi+p 2 + n + m\ p x ) \ p 2 



ACKNOWLEDGEMENTS 



(31) 



Research at Perimeter Institute is supported by the Government of Canada through Industry Canada and by the 
Province of Ontario through the Ministry of Research and Innovation. 



[1] P. Di Francesco, "2D quantum gravity, matrix models and graph combinatorics," math-ph/0406013. 

[2] P. Di Francesco, P. H. Ginsparg and J. Zinn- Justin, "2-D Gravity and random matrices," Phys. Rept. 254, 1 (1995) 
[arXiv:hep-th/9306153]. 

[3] B. Eynard, "Topological expansion for the 1-Hermitian matrix model correlation functions," JHEP 0411, 031 (2004) 
[hep-th/040726lj. 

[4] R. Gurau, "The 1/N expansion of colored tensor models," Annales Henri Poincare 12 (2011) 829-847, arXiv:1011.2726 
[gr-qc]. 

R. Gurau and V. Rivasseau, "The 1/N expansion of colored tensor models in arbitrary dimension," Europhys.Lett. 95 

(2011) 50004, arXiv:1101.4182 [gr-qc]. 

R. Gurau, "The complete 1/N expansion of colored tensor models in arbitrary dimension," Annales Henri Poincare 13 

(2012) 399-423, arXiv:1102.5759 [gr-qc]. 

[5] R. Gurau, "A generalization of the Virasoro algebra to arbitrary dimensions," Nucl. Phys. B 852, 592 (2011) 
[arXiv:1105.6072 [hep-th]]. 

[6] V. Bonzom, R. Gurau, A. Riello and V. Rivasseau, "Critical behavior of colored tensor models in the large N limit," Nucl. 

Phys. B853, 174-195 (2011). [arXiv:1105.3122 [hep-th] ] 
[7] V. Bonzom, R. Gurau and V. Rivasseau, "Random tensor models in the large N limit: Uncoloring the colored tensor 

models," Phys. Rev. D 85, 084037 (2012) [arXiv: 1202.3637 [hep-th]]. 
[8] R. Gurau, "Universality for Random Tensors," arXiv:1111.0519 [math. PR]. 

[9] M. Cho, D.Kim, S. Seo and H. Shin, "Colored Priifer codes for k-edge colored trees," The Electronic Journal of Combina- 
torics 11 (2004). 

[10] R.L. Graham, D.E. Knuth and O. Patashnik, "Concrete mathematics: a foundation for computer science," 2nd edition, 

Addison- Wesley, 1994. 
[11] R.P. Stanley, "Enumerative combinatorics, Vol. 2," Cambridge Univ Pr., 2001. 



